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Abstract. We simulate CMB maps including non-Gaussianity arising from cubic 
order perturbations of the primordial gravitational potential, characterized by the 
non-linearity parameter (/nl- The maps are used to study the characteristic nature 
of the resulting non-Gaussian temperature fluctuations. We measure the genus and 
investigate how it deviates from Gaussian shape as a function of (7nl and smoothing 
scale. We find that the deviation of the non-Gaussian genus curve from the Gaussian 
one has an antisymmetric, sine function like shape, implying more hot and more cold 
spots for > and less of both for ^nl < 0. The deviation increases linearly with 
(7nl and also exhibits mild increase as the smoothing scale increases. We further study 
other statistics derived from the genus, namely, the number of hot spots, the number 
of cold spots, combined number of hot and cold spots and the slope of the genus curve 
at mean temperature fluctuation. We find that these observables carry signatures of 
(7nl that are clearly distinct from the quadratic order perturbations, encoded in the 
parameter /nl ■ Hence they can be very useful tools for distinguishing not only between 
non-Gaussian temperature fluctuations and Gaussian ones but also between gNL and 
/nl type non-Gaussianities. 
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1. Introduction 

Since tlie discovery of the temperature anisotropies of the cosmic microwave background 
radiation [1] their statistical nature has been subject of intense study. If inflation 
[2, 3, 4], as strongly supported by observations, is indeed the mechanism that gave rise 
to these anisotropies, then their statistical nature must be inherited from those of the 
primordial density fluctuations. All models of inflation, in general, predict some amount 
of deviation of these fluctuations from a Gaussian distribution. The detailed knowledge 
of the deviations are quite model dependent. This makes it a good discriminant between 
various models of inflation. The observational search for non-Gaussianity is, however, 
beset with serious difficulties since various spurious observational effects can mask the 
true signal. 

Primordial non-Gaussianity arises from higher order terms in the perturbative 
expansion of the primordial gravitational potential, $, which must be taken into account 
in the presence of higher order interaction. In this paper, we consider the following 
expansion of the /oca/ type [5, 6, 7]: 

<l>(x) = <|.^(x) + /nl (($^(x))2 - ((<l>^)2)) + ^7nl($^(x))3 + _ _ _ ^ 

where /nl and (7nl are parameters that measure the level of non-linearity and <I>^ is the 
linear order perturbation. This expansion is rather special in that it assumes that the 
higher order perturbations are known in terms of the linear field, with our ignorance 
pushed into /nl and g'NL- Also, they depend only on the linear field value at the same 
spatial point. Such an expansion is not a generic prediction of all inflation models. In 
general, it can be a much more complicated expression involving convolutions of the 
products of "I'^'s {non-local) with the non-linearity parameters being scale dependent 
kernels. Generally, the predictions of inflation are quantified by n-point functions in 
Fourier space. Eq.(l) holds provided the n-point functions have most of the signal 
coming from some special configurations of the wave vectors. For the 3-point function, 
this corresponds to having the amplitude of one of the wave vectors tending to zero, the 
so called squeezed configuration. For the 4-point function, it corresponds to either one 
of the wave vectors tending to zero or two of them tending to zero [8, 9, 10, 11]. Here 
we ignore any possible scale dependence of /nl and (^nl to simplify the problem. This is 
justified by current experimental limitations. The bulk of the study on the topic, both 
theoretical and observational, has been focused primarily on the quadratic order non- 
Gaussianity [6, 12, 13, 14, 15, 16]. The strongest limits on /nl from CMB observations 
obtained so far is -4 < /nl < 80 (at 95% CL) [17]. 

There is now growing interest in the cubic order perturbation and its observability. 
The question is whether the cubic term is negligible, comparable or even dominant 
when compared to the quadratic term. For the standard single field and uncoupled 
mutiple fields slow-roll inflation they are found to be slow roll suppressed and hence 
negligibly small [9, 18]. (See [19] for attempt to get large trispectrum in single field 
inflation). However, in the curvaton scenario [20, 21] it can happen that the cubic term is 
comparable or even dominant to the quadratic term [22, 23, 24, 25, 26, 27]. Another class 
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of models with similar predictions is the so called multibrid inflation [28, 29, 30]. Thus, 
the possibility arises that the dominant source for primordial non-Gaussianity comes 
from the cubic non-linearity and it becomes important to understand theoretically how 
observational quantities in the CMB as well as the large scale structure get affected by it. 
Early studies of the consequences on the CMB have focused on the angular trispectrum 
of the temperature fluctuations [5, 31, 32], while feasibility study of measurement of 
the trispectrum from observational data was performed in [33]. Implications of (7nl on 
galaxy bispectrum was studied in [34]. The first limit on (^nl using SDSS data and 
N-body simulations has been obtained to be —3.5 x 10^ < g^i^ < +8.2 x 10^ (at 95% 
CL) in [35]. No bounds have been obtained using CMB data as yet. 

An important step towards understanding how the primordial non-Gaussianities 
affect the CMB temperature fluctuations is to simulate maps with the non-Gaussianities 
going in as input in the map making process. Since we have control over the input 
parameters, namely, /nl and (?nl in this case, we can test and calibrate the sensitivities 
of different statistical observables to these parameters, using the simulated maps. 
Observational contaminants can be added to find out exactly how each one of them 
can mask the real non-Gaussian effects and experimental bounds can be obtained for 
the input parameters by comparing with observational data. There are several map 
making methods that have been proposed in the literature for /nl non-Gaussianity. 
Komatsu et al. [36] used a straightforward method of generating a Gaussian random 
realization of the linear gravitational potential in Fourier space and then convolving two 
such fields to obtain the quadratic term in Eq.(lO). Liguori et al. [37] posposed a fast 
algorithm which require the computation to be done in real space rather than Fourier 
space. Other methods involve the input of given power spectrum, bispectrum and higher 
order spectra [38, 39], or some known correlation structure of the non-Gaussian field [40], 
or using sherical wavelets [41]. 

Simulated non-Gaussian CMB maps arising from (yfNL term have not been discussed 
in the literature as yet and we present them in this paper. The goal is to study how the 
cubic perturbations show up as non-Gaussianity of the temperature anisotropics. We 
use the genus, which is the number of isolated hot spots minus the number of isolated 
cold spots, as our statistical observable. We first make simulations of non-Gaussian 
maps with g-^i^ as the input parameter, extending the method of [37] and using the full 
linear radiation transfer function. In order to be able to investigate the pure effects 
of the (7nl term we have set /nl = in our simulations. Then, we use these maps 
to compute the genus statistic to find out how it varies as a function of (^nl and the 
smoothing scale. Further, we discuss four new statistics derived from the genus and 
show that they can be very useful tools to distinguish primordial non-Gaussianity from 
Gaussianity and also to distinguish between /nl and (^nl type non-Gaussianities. 

This paper is organized as follows: in section 2 we outline the method for generating 
the non-Gaussian maps, describe the implementation of the map making process and 
we present our results of the non-Gaussian maps and the one-point PDF. In section 
3 we compute the genus using the simulated maps and show how they deviate from 
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the Gaussian shape. We then discuss the derived statistics and elaborate on their 
characteristics and how they distinguish g-^i^ from /nl- We conclude with a summary 
of results and remarks on direction for future work in Section 4. 



2. Simulation of non-Gaussian maps 

We briefly review the method for simulating non-Gaussian maps outlined in [37] with a 
simple extension to include (?nl term. 



2.1. Calculating aim 's in real space 

The CMB temperature fluctuations are usually expanded in terms of spherical harmonics 
as AT{n) = J2im ozm^«m('T-)- The a^m's are then computed by convolving the primordial 
potential fluctuations with the radiation transfer function A^(r) as, 

r d^k 

a,m = M-^y J m Mk) Y;jk), (2) 

where $(k) is the Fourier transform of the real space potential $(x). Ai^k) encodes 
the evolution history of the CMB photons in their journey from recombination till now. 
Defining 

^im{k) = Jdnj^^k)Yem{k) . (3) 

we can rewrite a£m as 

(-iY C 

aim = / dk k"" ^im{k) Ai{k) , (4) 

To rewrite Eq. (4) as an integral in real space, define the real space harmonic 
potential 

^im{r) = ^ J dkk^^im{k)ji{kr) , (5) 
and its inverse 

^t^{k) = AniiY [ drr^^im{r)Ukr) , (6) 



where j^'s are spherical Bessel functions. Then, insert Eqn. (6) in (4) and define 

Aeir) = ^ Jdkk^Ai{k)ji{kr) , (7) 
we can then write: 



aem = J drr^^em{r)Ae{r) . (8) 

For non- Gaussian $ given by Eq.(lO), we would have 

<l>(k) = <|.L(k) + /NL$''^(k) + gNL'^^'^'H^), (9) 
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where 

d^k 



*^^'(k) = l^s§^s + ki + k.)$^(kO $^(k,). (10) 

We can then define harmonic components $^ and ^f^^ in Fourier and real space using 
Eqs. (5) and (6), to give us 

^,Uk) = ^Lik) + fNL^Zik) + gNL^Z'-ik). (11) 
aem is then given by 

O-em = (^Irn + /nL a^rn + 5'NL Ofo,^^: (12) 

where each term is an integral over the corresponding $. 

Thus, we need to compute four quantities, namely, Ae{r), $^^(r), $^(r) and 
$^^(r), in order to get aem upto cubic order primordial perturbations. Ai^r) can 
be independently computed using A£{k) obtained from CMBFAST [42]. In order to 
generate ^irni''^) need its correlation function given by [37], 

{<^imM)'^'^:mM)) = \^t^Z\ I dkeP^{k)jeAkr,)je,Akr2) , (13) 
where P<s>{k) is the primordial power spectrum (for the Gaussian part of $), given by 



with Ao being the amplitude, ko is some suitable pivot scale and is the spectral index. 
Then, ^frni''^) can be obtained from the integral 

^Ll'") = Jdrirlnem{ri)Wi{r,ri) , (15) 

where nim{f) are independent complex Gaussian variables characterized by the 
correlation function 

and We,{r^ri) are filter functions defined as 

W,{r, ri) = ^j dk e y^Mk)UkrMkr^) . (17) 

For fixed r, We{r, ri) is a smooth function of ri and sharply peaked at r = ri. Simplified 
expressions which are convenient for numerical calculation of Wi{r,ri) are given in 
Appendix A. 

To compute first compute the hnear potential $^(r) = J2em '^£m('^)^m('^) ^^^l 
square it to obtain $'^'"(r). Then harmonic transform to get $^(r). Similarly, 
can be computed by first taking cube of $^(r) to obtain <|)^'^^(r) and then harmonic 
transforming to get $^^^(r). Finally, putting /nl = we get, 

AT = AT^ + ^nlAT^^. (18) 

This method is particularly useful for calculating the (^nl case because had we done the 
calculation in k space we would have had to convolve three ^^^s as in Eqn.(lO). 
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Figure 1. Pixel Distribution of AT^'^ with respect to AT'^. 



2.2. Implementation of the algorithm 

We use a ACDM cosmological model with primordial spectral index = 1. We have 
used the WMAP 5 year parameters [43] given by VL^ = 0.233, = 0.0462, = 
0.721, = 0.087, /lo = 0.719. With these parameters, the conformal time today is 
To = 14360 Mpc~^. The accuracy of ^frni''^) t)e tested by computing its 'angular 
power spectrum' [44], which for = 1 is obtained as I'^Lil^)!^ ^ 1). 

We have sampled r at 472 points, with different step sizes chosen at different epochs, 
based on the shape of Ai{r). The accuracy of the resulting Gaussian Ce was tested 
by averaging over several maps and then comparing with the theoretical output of 
CMBFAST. We have used the Healpix package [45] to perform the harmonic transforms 
and the CMBFAST package [42] to compute the radiation transfer functions Ai[k). 
We have fixed £max = 1100 and the harmonic transforms have been computed using 
TVgide = 512. We find that the ^fNL term begins to dominate the linear term roughly 
around (7nl ~ 10^ and hence the perturbation expansion of $ is invalid beyond this 
value. The Gaussian maps are normalized by CMBFAST, while the non-Gaussian ones 
are normalized by matching the values of Ci at the first acoustic peak, i = 220, with 
the Gaussian one. 

In Fig.(l) we show for one simulation, how for a given pixel (^nlAT^*^ and AT*-' are 
correlated. Each of the green and blue dots represents a pixel. The left panel shows the 
pixel distribution of pure non-Gaussian temperature AT'^*~' about the (AT^)'^ curve. 
The right panel shows how the full AT, deviates from the AT^ line, for positive and 
negative values of (^nl- 
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Figure 2. All maps in this figure are obtained from one Gaussian realization. The left- 
hand ones are smoothed by FWHM = 7° while the right hand ones are smoothed by 
FWHM = 30'. The top ones are Gaussian, middle are non-Gaussian with gNL — 5x10^ 
and the bottom ones have ^nl = —5 x 10^. 

2.3. Non-Gaussian maps and 1-point PDF 

Gaussian and non-Gaussian maps obtained for the same Gaussian realization, for two 
different Gaussian smoothing scales, are shown in Fig. (2). The positive (7nl niaps show 
hot spots that are relatively hotter than those in the Gaussian map. On the other 
hand, the maps with negative (7nl show relatively cooler hotspots. The large (7nl value 
±5 X 10^ is chosen to make the differences in the maps visible. Note that the maps vary 
from realization to realization due to statistical fluctuations and at such large value 
of s^NL the non-Gaussian term of the temperature fluctuation may dominate over the 
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Figure 3. Effect of ^nl on 1-point PDF. Black (solid) line correspond to Gaussian, 
while green (dashed) corresponds to non-Gaussian with ^nl > and red (dotted) 
corresponds to jnl < 0. Results are averaged over 150 Gaussian and non-Gaussian 
maps smoothed with FWHM=30'. The right panel shows the deviation of P(AT) of 
the non-Gaussian maps from the Gaussian ones. 



Gaussian part for some realizations. We have shown here a reahzation which is close 
to the average behavior indicated by the average 1-point PDF described below and for 
which the non-Gaussian term is still sub-dominant to the Gaussian term despite the 
large value of (^nl- 

The 1-point PDF's, P(AT), are plotted in Fig. (3) for ^nl = ±5 x 10^ and 
FWHM=30', averaged over 150 realizations. AP is the difference between non-Gaussian 
and Gaussian PDF's and we have shown them in the right hand side of Fig. (3) for 
positive and negative S'nl's. We see that for both positive and negative (^nl, the mean 
position is not affected, as expected from the fact that non-Gaussian part comes from 
($^) . Positive (7nl increases the pixels around the mean temperature fluctuations, 
decreases the intermediate temperature range and again increases the hottest and coldest 
ranges, leading to leptokurtic shape of the 1-point PDF relative to the Gaussian one. 
Negative (7nl has the opposite effect and results in platykurtic shape of the 1-point PDF 
relative to the Gaussian one. These effects become more pronounced as we increase 
(^NL and corroborates what we observe visually in the maps shown in Fig. (2). On the 
angular power spectrum, C^, both positive and negative g^i^ increase the low scale power 
with negative (^nl having a stronger effect. 

3. Genus statistic 

By means of iso-temperature contours of the temperature fluctuation field one can study 
its global morphological properties. The genus, which is the number of isolated hot spots 
minus the number of isolated cold spots, can be obtained from the contours for a given 
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Figure 4. Variation of g(y) with ^nl for different smoothing angles . Results have 
been averaged over 200 realizations. G<3,max .(.j^^, amplitude of the Gaussian genus 
evaluated aX,v = \. Black (solid) lines are Gaussian, green (dashed) lines correspond to 
postive 5nl and red (dotted) to negative 5nl- The values of 5nl are ±5 x 10^, ±1 x 10^ 
and ±2 x 10^. The shaded regions in the top panels show the l-cr error bars for the 
Gaussian genus curve. 



threshold temperature, denoted by z/ = AT/fToj where ctq is the standard deviation of 
the temperature fluctuation. It is sensitive to the Gaussian/non-Gaussian nature of 
the fluctuation field and the shape of the underlying angular power spectrum. This 
makes it a useful tool to test non-Gaussianity. It was introduced in the context of 
the CMB in the seminal papers [46, 47]. It has been used extensively to study non- 
Gaussianity in a number of papers [48, 49, 50]. The genus is one of the three Minkowski 
functionals (MF's) [47, 51, 52], which are topological quantities that can be defined 
for a two dimensional field and which completely characterize its topological properties. 
The other two are the total iso-temperature contour length and the fraction of total area 
above the threshold. They have been used to constrain the /nl parameter [43, 53, 54, 55]. 
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For a given temperature threshold the genus is given by 

G{v) = — [ nds, (19) 



'271 jc 

where k is the signed curvature of the iso-temperature contours G. The genus can also be 
parametrized by the area fraction above the threshold. Using the temperature threshold 
is computationally easier, while the area fraction decreases correlations between the 
MF's [56]. Here we use the temperature threshold since we are focussing on the genus 
only. 

For a given Ge and Gaussian smoothing angle 6s, related to FWHM as 6s = 
FWHM/2v21n2, the genus per steradian of a Gaussian temperature field can be 
expressed as 



1 E^(^ + i)(2^+i)c^^^ 



2 



where Fi is the smoothing filter, which for a Gaussian filter is given by Fi = exp(— + 
l)6'^/2). For a weakly non-Gaussian field characterized by /nl, approximate analytic 
expressions for the Minkowski functionals were obtained in [57, 58]. The expressions 
upto (7nl order on the perturbations have not been calculated yet. Here we study the 
topology of the temperature fluctuations that arise purely from (jfNL non-Gaussianity by 
measuring the genus. For the computation of the genus we follow the method of [47]. 
We have used 31 threshold levels in the range -3 to 3. The amplitude and shape of 
the genus curve is sensitive to the smoothing scale and we have chosen three scales, 
FWHM = 30', 45' and 60' to demonstrate our results. Results become more noisy at 
higher smoothing scales since the number of structures (hot and cold spots) decrease. 
We have averaged the results over 200 maps for each value of (7nl- 
Let us denote: 

AG{iy)=G''^{u)-G^{u), (21) 

where G'-^^u) is the Gaussian genus and G^'^^u) is the non-Gaussian one. In Fig. (4) we 
have plotted G and AG/G"^'™'^^, where G'^''^^^ is the amplitude of the Gaussian genus at 
v = 1. We have shaded the region within l-cx error bars for the Gaussian genus curves 
in the upper panels. It shows that for positive (^nl, the amplitude of the non-Gaussian 
genus curve is higher than the Gaussian one in the threshold range < < 1, while 
it is lower in the range 1 < \v\ < 2.5. The fact that the genus is smaller in the range 
1 < |z/| < 2.5 means there are fewer hot spots and cold spots in the CMB map when 
non-Gaussian contribution with positive (^nl is present. Because of the larger genus 
amplitude at < < 1, the range of u showing the sponge-like topology [59] (actually 
a two dimensional cut through sponge in the context of CMB) is smaller relative to the 
Gaussian one. When g^i^ is negative, there are more hot as well as cold spots. Other 
than an overall scaling of the amplitude of AG, the smoothing scale does not seem to 
affect its shape in the threshold range 1 < \v\ ^2.5. 

In order to quantify the functional dependence of AG/G*^'™^^ we take the average 
of the magnitudes of its peak values (which lie roughly in the ranges 1.5 < < 2 and 
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Figure 5. Functional dependence of AG/G^'"^'^^ on ^nl and smoothing scale. In 
the right hand figure, the lines from the bottom to the top correspond to ^nl = 
±5 X 10^±106,±2 X 10^ . 



— 2 < 1/ < —1.5). Let us denote it by AG. In the left panel of Fig. (5) we have plotted 
the variation of AG with gfNL for FWHM 30', 45' and 60'. We find a linear dependence on 
qtssl, indicating that the leading order contribution to the deviation from the Gaussian 
genus curve comes from terms of order (^nl- The right panel shows the variation of AG 
with the smoothing scale for fixed (^nl- We find a mild increase as we increase FWHM, 
in the range of smoothing scales that we have studied. 

So far we have described the number of structures for threshold range |z/| < 2.5. 
Upto the smoothing smoothing scale of FWHM= 60' the number of structures above this 
threshold are not significantly large. However, an interesting observation that we have 
made for smoothing scales above FWHM= 60' is that the relative number of structures 
in this hottest or coldest range of threshold values, > 2.5, compared to the range 
1 1^1 < 2.5, grows significantly and hence promises to be useful for constraining (/nl at 
higher smoothing scales. We have not shown the results since the plots are quite noisy 
but we will be exploring this region further in subsequent work. 

3.1. Other statistics derived from the genus 

The genus at different threshold values is strongly (anti) correlated. One may then think 
of deriving other observables using the information inherent in the genus curves so that 
the non-Gaussian information may be maximized, and which may distinguish between 
different kinds of non-Gaussianity. The simulated non-Gaussian maps may then be 
used to test the sensitivity of these observables. We mention here four such quantities, 
namely, Rcoid, R-hot, Rspots and So, which are defined below. Then we explore how they 
deviate from the Gaussian expectations. 
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Figure 6. Rhot vs. Rcoid for different values of ^nl for FWHM = 30'. The open circle 
at the center is the Gaussian mean while filled ones are ^nl = ±5 x 10^, ±1 x 10^ and 
±2 X 10^, at increasing distances from (1,1). The region where the mean values for 
positive and negative values of /nl must lie, provided FWHM < 94', are indicated. 
For FWHM > 94', they will interchange quadrants. 



(i) Rcom: Let Ncom denote the total number of cold spots, defined as Ncom = 

du G{u), where i^i, 1^2 are suitably chosen positive threshold values with U2> ui. 
Let G^^iv) denote the Gaussian curve obtained by fitting the non-Gaussian genus 
points at different threshold values to a Gaussian shape. Let N^j^ = duG^^{u). 
Then, we define 

^^cold 

For Gaussian maps, its value must be one. By inspecting Figs. (4) and choosing 
Pi = 1 and U2 = 2.5, we can predict that for (^nl > 0, Ncom must be less than one, 
whereas, for gjssL < it must be greater than one. 

(ii) Rhot: is defined to be the ratio Nhot/N^^^, where Nhot and ^^^t defined similar 
to Ncoid and N^j^, but for integration in the positive threshold range. Again the 
Gaussian expectation is one, and for same ui and 1^2 as above Rhot < 1 for (^nl > 
and Rhot > 1 for (^nl < 0. If AGi^u) is antisjnumetric, as it appears to be in Fig. (4), 
then we must have Rhot = Rcoid- 

(iii) Rspots: is defined to be the ratio Ngpots/N^^ts' where N^pots = N^oid + Nhot and 
N^ots = ^mid + Nhot- The Gaussian prediction of this statistic is one. For ui = \ 
and V2 = 2.5 it is less than one for (^nl > and greater than one for (^nl < 0. 

(iv) So : is defined as the ratio of the slope of the non-Gaussian genus curve and that 
of the fitted Gaussian at = 0. It is greater than one for (^nl > and less than 
one for g^i^ < 0. 
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In the case of non-Gaussianity arising from /nl, since the shape of AG is strongly 
dependent on the smoothing angle, these observables will depend on the smoothing 
angle. By judicious choice of vi and V2 "we can make them maximize the differences 
between /nl and (7nl type non-Gaussianities. By inspecting Fig. (2) of Ref. [58] and 
choosing Ui = \ and z/2 = 2.5, one can deduce that we must have Rcow < 1, R-hot > 1 
for positive /nl and Rcoid > 1, Rhot < 1 for negative /nl, when FWHM < 94' . The 
situation is reversed if FWHM > 94'. Sq is always greater than one for positive /nl and 
always less than one for negative case. For FWHM < 94', Rspots > 1 if /nl > and 
Rspots < 1 if /nl < and they interchange for FWHM > 94' . 

We next analyze the properties of these four observables inferred from the non- 
Gaussian maps. An overhead bar denotes the mean value of each of the above four 
statistics obtained by averaging over 200 realizations. Fig. (6) shows the parameter 
space (Rhot, Rcoid) for varying (^nl values, for FWHM= 30'. We have used z/i = 1 and 
V2 = 2.5 to calculate them. The open blue triangle at (1, 1) indicates the Gaussian 
mean, while the red filled ones denote (^nl = ±5 x 10^, ±1 x 10^ and ±2 x 10®, at 
increasing distances from (1,1). In practice, the Gaussian mean has a small shift from 
(1,1) due to statistical fluctuation. We have corrected it by simply shifting it. We have 
then shifted the non-gaussian means by the same amount. This shift does not affect the 
relative distance between Gaussian and non-Gaussian means. Note that with minimal 
statistical fluctuation and if AG is exactly anti-symmetric then the slope of the line 
(Rhot, Rcoid) will be exactly 45°. 
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Figure 8. Functional dependence of Rspots on jinl and smoothing scale. In the right 
hand figure, the lines get farther away from Rspots — 1 Ihie as I^nlI increases. The 
values are c^nl = ±5 x 10^, ±10*', ±2 x 10^. 



The parameter space (Rhot, Rcom) gets divided into four regions, with negative ^^nl 
occupying the first quadrant and positive ^^nl the third, while =f/nl occupies the second 
and fourth quadrants, respectively. Larger smoothing angle increases the dispersion of 
the distribution of individual points (Rhot? Rcom) about the mean because the number 
of structures decreases. But at the same time the size of deviations increases as 
smoothing scale increases, making the low resolution maps also useful in discriminating 
non-Gaussian maps from Gaussian ones. 

Next we discuss the remaining two observables, Rspots and Sq. Fig. (7) shows 
the parameter space spanned by them for varying (yfNL- and V2 are the same as 
above. Again the blue triangle indicates the Gaussian mean, while the red ones denote 
g^i^ = ±5 X 10^, ±1 X 10^ and ±2 x 10^, at increasing distances from (1,1). The 
smoothing scale is FWHM = 30', and the effect of smoothing scale is similar to above. 
The parameter space gets divided into four regions, with if^NL occupying the second 
and fourth quadrants, respectively, while ±/nl occupies the first and third quadrants, 
respectively. Again for FWHM > 94', the positions of +/nl and — /nl will interchange. 

Since Rspots carries information of both Rhot and Redd we choose it for showing 
functional dependence of (7nl and FWHM. Panel 1 of Fig. (8) shows how Rspots varies 
with (7nl when FWHM is fixed. As is clear from the figure we find linear dependence, 
which agrees with the fact that we saw linear dependence of AG on g^ssL in Fig. (5). 
Panel 2 of Fig. (8) shows functional dependence on FWHM for fixed ^^nl- The lines 
get farther away from the Gaussian expectation, Rgpots = 1, as we increase |(7nl|- They 
exhibit mild increase with increase of the smoothing scale. 
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4. Conclusion 

We have simulated non-Gaussian CMB maps with the non-Gaussianity coming from 
purely third order perturbations of the primordial gravitational potential. We have 
used the map making algorithm proposed by [37] which computes a^m's as an integral 
in real space. This method is particularly advantageous for including the third order 
linearity, as compared to integrating in Fourier space since there are two k convolutions 
involved which make the computational very heavy. 

We have investigated how the 1-point PDF gets modified from Gaussian shape due 
to the effect of the non-linearity parameter (^nl- We found that positive (?nl changes 
the 1-point PDF to leptokurtic shape and negative (^nl changes it to platykurtic shape. 
Its effect on the C^'s is to increase power for large £'s for both positive and negative (^nl 
with the effect of negative (7nl being stronger. 

We have next used the simulated maps to compute the genus curve and calculate 
their deviation from the Gaussian shape. The purpose is to understand their functional 
dependence on g^i^ and how the non-Gaussian term modifies the topology of the CMB 
temperature field. We found that positive (/nl decreases both hot and cold spots in the 
threshold range 1 ^ ^ 2.5 and increases the genus in the range < |z/| < 1. The 
effect of negative (^nl is just the opposite. This results in antisymmetric shapes of AG, 
which look approximately like sine functions. We do not find significant variation of the 
overall shape of AG as we vary smoothing scale, but the amplitude of AG decreases 
as we increase the smoothing scale. We found that the sensitivity of the negative and 
positive (^nl's are roughly same. These results are very different from the genus arising 
from /nl, which leads to symmetric form of AG and which has strong dependence on 
the smoothing angle. 

We have also studied four other statistics derived from the genus, namely, the 
number of hot spots, the number of cold spots, the combined number of hot and cold 
spots and the slope of the genus curve at z/ = 1. We found that these quantities 
carry distinct signatures of /nl and (?nl- The parameter spaces get neatly divided into 
quadrants with each of the positive and negative (^nl and positive and negative /nl 
occupying one. Hence they can be very useful for distinguishing these two different 
types of non-Gaussianities. 

Since our goal was to get theoretical understanding of the nature of non-Gaussianity 
arising from (yfNL term by means of studying statistical observables such as the genus and 
other quantities derived from it, we have not considered real observational contaminants 
such as point sources, instrument noise etc., to our simulations. One needs to take them 
into account for actual comparison with experimental data and putting constraints ^^nl- 
Another interesting observable that we are studying using the non-Gaussian simulations 
is the correlation of peaks in the maps. These will be the subject of forthcoming 
publications. 
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Appendix A. Simplified expressions for Wi 

Here we present expressions which simplify the numerical computation of the filter 
function given in Eq. (17). Instead of doing the integral over product of two j^'s, which 
is very time consuming due to the highly oscillatory behavior of j^, we can simplify 
Wi to express it in terms of Gamma functions. Using je{kr) = ^J 7r/2krj£^i/2{kr) and 
P$(/c) =Ao^we get, 

Wi{r,r,) = ^J^J^°^ dkk-'^Me+,/,{kr)Je+,/,{kr,) (A.l) 

CASE 1 : r = ri. Using the following formula [60], 

dkk-^J,{ak)J^{ak)- y 2 ) 



X ^nz^, (A.2) 

which holds when i/ + /i + l>A>0 and a > 0, and putting A = 1 — ns/2 we get, 

w( ^- ^ 1 r(i-f) + \ + 

yy^r. rij ^1-^ ^1+^ [r(i _ i^)]2 Y{i + i-^y ^""-^^ 

CASE 2: r ^ r\. For this case, when = /x, we can use the following formula [60], 
dk k-^J,{ak)J,{f]k) - ' y ^ 2 ) 



2^(a + /3)2--^+i r(z/+ i)r(i±^) 



which is holds provided 2z/ + 1 > A > — 1 and a > 0, /? > 0. -F(a, 6; c; z) is the Gauss 
Hypergeometric function given by 

Tie') fb-i(i _f]c-b-i 

^f"- ^' ^) = fWi^) I * (1-4 • 

which is valid when c > 6 > 0. Then, using 

T{2i) = -1=22^-1/2 p(^) r + (A.6) 
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and defining z as 

_ 4rri 4ri/r 
(r + riY (1 + ri/rY ' 



we get 



' V TT (r + ri)iW2 r(l-i^) r(£+l 



2 ' 4 



where 



(A.7) 



X / c^t — ^ -r^- (A.8) 



We can then scale the r dependence as, 

^^(^' ^i) = ;:t^ ^^(n/^), (A.9) 

WArJr) = 2"^/' - r(f + 1 + ^) 

^^1/^ Vtt (l + ri/r)i+W2 r(l-f) r(£ + l) 

(1 — tZ)^^2+ 4 

Since r 7^ ri we have < 2; < 1. The integrand in Eq. (A. 10) is a smooth positive 
function with a local maxima in the interval [0:1]. It becomes more and more localised 
as i increases and the peak position shifts towards one as z approaches one. The integral 
can be easily computed numerically. 
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